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ABSTRACT: Pulsed laser deposition is one of the most promising modern approaches for producing 
epitaxial layers. The approach gives a possibility to apply materials with special properties (metals, 
carbides, etc.) to the surface of parts, which allows restoring geometry, increasing surface strength and 
corrosion resistance, etc. This paper considers an analytical approach for solution of partial differential 
equations. This approach has been used for analysis of mass and heat transfer in reaction chamber during 
growth of the epitaxial layers by using pulse laser deposition. Changing of mass and heat transfer is 
investigated depending on a number of parameters. We consider an analytical approach for analysis of 
mass and heat transfer. The approach gives a possibility to take into account changing of parameters of 
processes both in space and time, as well as the nonlinearity of the processes. The approach was used for 
the analysis of mass and heat transfer in the reaction chamber during growth of epitaxial layer by using 
pulsed laser deposition. Changing of mass and heat transfer was investigated depending on series of 
parameters. The approach for analysis mass and heat transfer, which were considered in this work, makes 
it possible to carry out a more adequate prognosis of growth of an epitaxial layer by using pulsed laser 
deposition in comparison with similar approaches. 
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INTRODUCTION 
One of the most promising methods of growing : OT (x,t) (Tu) _e a(n)? T (x,t) ey (1) 
of epitaxial layers is a method based on pulsed lar ax | ax ax |?” 


laser deposition. This method gives a possibility to 
grow layers of materials with special properties 
on the surface of details. The application gives a 
possibility to restore geometry, increase surface 
strength and corrosion resistance, to manufacture 
structures with the required properties [1-10]. In 
this paper, we consider mass and heat transfer in a 
reaction chamber during growth of an epitaxial 
layer by pulsed laser deposition. We introduced an 


where gis the density of the evaporated material; 
Cp is the specific heat at constant pressure; T (x,t) 
is the temperature distribution; 2 (T) is the 
thermal conductivity (dependence of thermal 
conductivity on temperature could’ be 
approximated by the following relation: / (T) 
=Aass{1t+u [Ta/T(x,t)]%}, see, for example, [11]); p 
(x,t) is the power density of laser radiation; x and t 


analytical approach for predicting of the are the current coordinate and time; a (T)=A 
considerate process. The approach gives a (T)/c (T) is the thermal diffusivity. Speed of 
possibility to take into account nonlinearity of the movement of the evaporation boundary is 
considerate processes as well as changes of determined by flows Ji of particles evaporated 
parameters of the process in space and time. from the surface: U (t)=>J,/p,, where i 


means the material used during growth. Boundary 
and initial conditions could be written in the 
following form 


T (:0,t) =T,, T (x,0) =Tr. 


= QO, “Uu (t)’ 


0) x 
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Fig-1 The direction of motion of the material Here 7; is the equilibrium temperature equals 


vaporized during laser deposition 
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room temperature; Q, is the heat of vaporization. We describe the transfer of the growth components using 


the second Fick law in the following form 


ACU) _yFClus) _ ae 0 C(x,t) - 
Ot Ox Ox Ox 
with boundary and initial conditions 
C (0,t) =Co, C (,t) =0, C (0,0) =Co, C (x >0,0) =0, (2a) 


where C (x,t) is the concentration of vaporized 
material; Dc is the diffusion coefficient of 
considered material. The concentration 
dependence of the diffusion coefficient could be 
written as D =Dass{1+é& [C (x,t)/P(x,T)]} (see, for 
example, [12]), where P(x,T) is the solubility of the 


corrections [13]. Framework this method, we 
replace the unknown functions T (x,t) and C (x,t) 
by their unknown average values air and aic in 
the right-hand sides of the considered equations. 
Then we obtain the equations for the first-order 
approximations of the considered functions 71(x,t) 


vaporized material. Next, let us to solve Eqs. (1) and C: (x,t) 
and (2) by method of averaging of function 
aT(xt)_¢ OC (ast) _ 
; = p(3,t), = (3) 
Ot Ot 
Solutions of the above equations could be written as 
t 
c,pT(x,t)=J p(x%,7)dt+T.,C(~1)=C,. (4) 
0 


The second-order approximations of the 
required functions T (x,t) and C (x,t) have been 
determined by the framework of standard 
procedure [13,14], i.e. by replacing these functions 
by the following sums: T (x,t)>aeartT1(x,t) u C 
(x,t) aect+Ci(x,t) on the right side of equations 
(2), where azr and azc are the average values of 


considered temperature T2(x,t) and concentration 
C2(x,t). Higher-order approximations could be 
calculated similarly with a corresponding 
increasing of number of order of approximation. 
The relations for the second-order approximations 
of the considered functions after the above 
substitution take the following form 


the second-order approximations of the 
at(ut)_, ([,, at lariat) 
. Ot ee & [a,,+T (x,t)? | Ax 
+c,pu (jhe!) + p(x,t), (5a) 
age 
AC,(x,t) _ a py eect Glan 2 Cut) u(t) C eet) (5b) 
Ot OX P’(x,T) Ox Ox 


Integration of the left and right sides of equation 


(5a) and (5b) on time gives a possibility to obtain the 


following relations for required approximations of the considered functions 


4) 2th, ate aren) 
epine)=2, =| ater P os dtt+ 
+c,pJu (72h) 4 r+{p(x,r)dr+c,pT,, (6a) 
0 Ox 0 
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OC, (x,r) 


Ox x 


Average values of the second-order approximations of considered functions a@zrand azc could be calculated 


by using the following standard relations [13] 


] oL 1] eL 
sp = Bn) -Tenld xde, t= Ble )-CGkd)ld xa. (7) 


Substitution of relations (4) and (6) into relations (7) gives a possibility to obtain the following results 


1 @ uf A O UT? 6 T,(x,t) 
—___f(@- ‘ass 1 d 1 
ea az! z C, pox "Te, +h) Ox * 
+c, pu (yt) dxdt, (8a) 
Cx 
1 @ L O [a,. +C,(x,t)f 2 C,(x,t) 
——__ /((@~- D —2]1 2C 1 1 
=a yal zl oO I a6 P’(x,T) Ox 
+u (Za) dxdt. (8b) 
Ox 
The average values of azr and azc depend on DISCUSSION 


values of parameters 42 andg. Framework this 
paper, we calculate concentration of the growth 
component and temperatures are as the second- 
order approximations framework the method of 
averaging functional corrections. The 
approximations are usually enough adequate 
approximations for obtaining qualitative 
consequences and obtaining some quantitative 
results. The results of analytical calculations were 
verified by comparison of them with results of 
numerical simulation. 


1.00 


0.00 L/4 


In this section, we have analyzed spatio- 
temporal distribution of concentration of 
vaporized material. Figs. 2 and 3 show 
dependences of concentration of vaporized 
material of frequency of laser pulses and 
temperature of substrate at fixed value of other 
parameter. The increasing of this concentration is 
natural, because with an increasing of the 
temperature and frequency of the laser pulses lead 
to increasing of quantity of evaporated material. 
Fig. 4 shows dependences of concentration of 
vaporized material on pressure in growth 
chamber. 
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Fig. 2. Distributions of the concentration of the growth component at various values of the laser pulse 
frequency. Increasing of curve number corresponds to increasing of the laser pulse frequency 
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Fig. 3. Distributions of the concentration of the growth component at various values of temperature of the 
substrate. Increasing of curve number corresponds to increasing of the temperature 
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Fig. 4. Distributions of the concentration of the growth component at various values of pressure in growth 
chamber. Increasing of curve number corresponds to increasing of the pressure 


CONCLUSION 


We consider an analytical approach for analysis 
of mass and heat transfer during pulse laser 
deposition in a reaction chamber during growth of 
epitaxial layer. The approach gives a possibility to 
take into account changing of parameters of 
processes both in space and time, as well as the 
nonlinearity of the processes (in this situation the 
approach could be consider as more adequate 
prognosis in comparison with previously 
considered approaches). We analyzed changing of 
mass and heat transfer in the reaction chamber 
under influence of variation of parameters. 
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